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ABSTRACT 

We report here on the longest deep X-ray observation of a supergiant fast X-ray transient 
(SFXT) outside outburst, with an average luminosity level of 10'^'^ erg s^^ (assuming 3 kpc 
distance). This observation was performed with Suzaku in December 2009 and was targeted 
on IGR J08408^503, with a net exposure with the X-ray imaging spectrometer (XIS, 0.4-10 
keV) and the hard X-ray detector (HXD, 15-100 keV) of 67.4 ks and 64.7 ks, respectively, 
spanning about three days. The source was caught in a low intensity state characterized by 
an initially average X-ray luminosity level of 4x10''^ erg s^^ (0.5-10 keV) during the first 
120 ks, followed by two long flares (about 45 ks each) peaking at a flux a factor of about 
3 higher than the initial pre-flare emission. Both XIS spectra (initial emission and the two 
subsequent long flares) can be fitted with a double component spectrum, with a soft thermal 
plasma model together with a power law, differently absorbed. The spectral characteristics 
suggest that the source is accreting matter even at this very low intensity level. From the HXD 
observation we place an upper limit of 6x10'^'' erg s^^ ( 15-40 keV; 3 kpc distance) to the hard 
X-ray emission, which is the most stringent constraint on the hard X-ray emission during a 
low intensity state in a SFXT, to date. The timescale observed for the two low intensity long 
flares is indicative of an orbital separation of the order of 10^^ cm in IGR J08408-4503. 
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1 INTRODUCTION days has been suggested based on the duration of the flaring activ- 
ity ( Romano et al. 2009 ). although a secure determination based on 

IGRJ08408~4503 is a hard X-ray transient discovered in 2006 ^ ,, .■ i j i .■ .n i i ■ ai i 

„ „ , ^ X-ray and/or optical modulation is still lacking. Also a pulse pe- 

May with /iVr£;G7?AL during a flare with a duration of 900 s, ..f u * j .u .u . -p.u . x.- . ■ 

. „ „ r- n nod has never been round, thus the nature oi the compact obiect is 

with a flux of 250 mCrab (at peak, 20-40 keV) JGotz et all , f j 

I2OO6I) . Archival INTEGRAL observations of this field before "" °^"' 

the discovery observation demonstra ted the recurrence of its flar- IGR 108408^503 has been classified as a member of the class 

ing activit y (Mereghetti rtalj \200di) . A refined X-ray position of high mass X-ray binaries (HMXBs) called supergiant fast X- 

JKennea & Cam pana 200^, allowed the association of the X-ray ray transients (SFXTs), X-ray transients with bright and short du- 

transient with a bright 08.5Ib(f) star ( HP 74194) located in the ration flaring activi t y associated with blue supergiant companions 

Vela region, at a distance of about 3 kpc JMasetti et a l] |200d) . Op- JSguera et aUl2005l Ismith et alj|2006l) . The large majority of the 

tical spectroscopy of HD 74194, performed just a few days after members of this class have been discovered by INTEGRAL 



the source discovery, showed variability in the Hq profile and a during the survey of the Galactic plane JSguera et afl 1200^ 

radial velocity variation with an amplitude of abo ut 35 kms"^ ISguera et al]|2006n . A few SFXTs are X-ray pulsars, thus demon 

(Hel and Hell absorption lines; iBarbaet al.ll2006l) . Other bright strating that the compact object is a neutron star, while in the other 

X-ray flares were observed with INTEGRA L and Swift in sources a black hole cannot be excluded, although the broad band 

October 2006 , July 2008 and S epteni ber 2008 JGotz et alj r 2007l . spectra (0.1-100 keV) in outburst are very sim ilar to those typi- 

iRomano et ai][2009l ISidoli et ij]|2009l) . An orbital period of ~35 cally observed in accreting pulsars (e.g. lRomano e t al. 2008 ). Their 

outbursts typically show a duration of a few days dRomano et al.l 

I2OOI ISidoli et al.ll2009h , with s niall duty cycles alth ough highly 

* E-mail: sidoli@iasf-milano.inaf.it variable from source to source JRomano et al J 120091) . Each out- 
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burst is characterized by several short (a few minutes to a few 
hour duration) and brig ht flares, each rea ching a luminosity of a 
few lO^^-lO^'^ erg s"^ JSidoli et alJ(2009 V A much lower inten- 
sity emission characterizes the long term beh aviour of these tran- 
sients (10^^-10^" erg s~ Msidoliet^l2008h . while the quiescent 
state (no accretion) has been observed rarely and displays a very 
soft (likely thermal) spectrum with a luminosity of '^lO'^^ erg s~^ 
jin't ZanduOOSllLev der et al. 2007). A very low level of accretion 
(~10^^ erg s-i) has been observed from IGR J08408-4503 in 2007 
with XMM — Newton as well, displaying a double component 
spectrum composed of a thermal soft part (likely from the super- 
iant companion), together with a harder power-law component 
Bozzo et alj2010t) . 

Although SFXTs only sporadically undergo a bright outburst 
(from null to 4 outbursts per source have been observed dur- 
ing the two y ears long monitoring of four SFXTs with Swift, 
iRomano et al.l 12009). their properties during outbursts are much 
better studied than their behaviour during the long-term lower in- 
tensity states (intermediate level or the quiescence). 

We report here on a Suzaku observation of IGRJ08408- 
4503 with a net exposure of 67.4 ks (XIS), spanning about three 
days. This is the longest X-ray observation of a SFXT catching the 
source in a very low intensity level. 

We also report on the analysis of our private INTEGRAL 
observations of the source field covering about one month including 
times of the Suzaku observation. We then compare our Suzaku 
results with a much shorter observation of a low intensity state from 
IGR J08408~4503 performed with XMM - Newton in 2007. To 
this aim, we re-analysed the se XMM — New ton observations, 
although already reported bv lBozzo et al.l J2010l) . 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Suzaku 

A Suzaku JMitsuda et ^ l2007h observation of IGRJ08408- 
4503 was performed from 2009 December 11, starting at 17:53 
UT, until December 14 12:00 UT. The target was placed at the 
'HXD nominal' pos ition. The X-ray Imaging Spectrometer (XIS; 
iKovama et alj2007h was operated in the normal mode with no win- 
dow option, p roviding a time resolut ion of 8 s. The Hard X-ray De- 
tector (HXD uTakahashi et al.ll2007r) was in standard mode, with a 
time resolution of 61 fis for the individual events. 

The XIS and HXD data were both processed with the Suzaku 
pipeline version 2.4.12.27 in the HEASOFT (version 6.6) software 
package. The data were screened according to standard criteria. In 
particular, events were discarded if they were acquired during pas- 
sages through the South Atlantic Anomaly or in regions of low ge- 
omagnetic cut-off rigidity «6 GV for XIS and <8 GV for HXD), 
or with low Earth elevation angles; for the XIS only events with 
grade 0, 2, 3, 4, and 6 were considered, and the CLEANSIS script 
was used to remove hot or flickering pixels. The net exposures ob- 
tained with the XIS and HXD (live time 93.2%) were 67.4 ks and 
64.7 ks, respectively. 

The XIS events of IGR J08408-4503 were accumulated in 
each of the three XIS cameras within a circular region (2' radius) 
centered on the target, while the backgrounds were estimated from 
source-free regions well outside the point-spread function of the 
source (see Fig. [T). For the spectral analysis, spectral redistribu- 
tion matrices and ancillary response files for the XIS spectra were 
generated using the XISRMFGEN and XISSIMARFGEN tasks. 
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Figure 1. XISl image of the IGR 108408^503 field. The two circular re- 
gion.s mark the extraction regions for the source and the background. In 
the upper pa rt of the X-ray image, di ffuse emission from the Vela SNR is 
evident (e.g.. lLu & AschenbachlboOOh . 



Background subtraction for the HXD data is performed us- 
ing a synthetic model that accounts for the time-variable particle 
background (the so-called 'non X-ray background', NXB). In par- 
ticular, we used the 'tuned' NXB model (see Fukazawa et al. 20()9| 
for details). The cosmic X-ray background and possible Galactic 
diffuse emission remain in the background subtracted data, as well 
as emission from any X-ray source that occurs in the instrument 
34 X 34 arcmin^ field of view. 

After subtracting the NXB, a positive signal is detected in 
the HXD-PIN up to ^40 keV (about 10% of the total counts in 
the 15^0 keV band); no significant emission is detected in the 
GSO data. To account for the expected contribution from the 'cos- 
mic X-ray backg round' (CXB ), we use d the spectra reported in 
iGruber et all ( [1999) and in Moretti et all J2009) . We estimate that 
the CXB contributes to the observed signal with roughly half of the 
NXB -subtracted counts. 

The spectrum of the remaining signal is adequately described 
(xl = 1.16 for 64 degrees of freedom) by a power-law model with 
photon index LOlpg (la confidence level) and 15^0 keV flux of 
~6 X lO"'^^ erg cm~'^ s~^. Considered that the combined uncer- 
tainties on the CXB and NXB are comparable to the flux estimated 
in this way, and lacking any specific signature of IGR J08408-4503 
(such as it would be a periodic signal) in the HXD-PIN band, we 
regard the observed flux only as an upper limit on the hard X-ray 
emission of the source. This upper limit is not very constraining 
for the extrapolation of the XIS spectrum at high energy. Thus, in 
the following, we concentrate on the soft energy part of the source 
spectrum observed with XIS. 



2.2 INTEGRAL 

The INTEGRAL (INTErnational Gamma-Ray Astrophysics Labo- 
ratory, Winkler et al. 2003) satellite, launched in October 2 002, car- 



ries 3 coded-mask telescopes : IBIS (15 keV- 1 Me V; Ubert ini et al 
|2003| ), the spectrometer SPI (20 keV-8 MeV; .Vedrenne et al 
200J), two X-r ay monitors JEM-XI and JEM-X2 (3-35 keV; 



Lundet al J 120030 . IBIS is composed of two detector layers: IS- 
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Table 1. Summary of the INTEGRAL data we have analysed. 



Rev. 



Period (MJD) 



Date 



863 55142.852-55144.414 

864 55146.289-55147.531 

865 55147.637-55148.500 

866 55150.715-55152.578 

867 55155.910-55156.500 

868 55156.613-55158.035 
880 55194.098-55195.400 



2009-11-07 20:26 - 2009-11-09 09:56 
2009-11-11 06:56 - 2009-11-12 12:44 
2009-11-12 15:17 - 2009-11-13 12:00 
2009-11-15 17:09 - 2009-11-17 13:52 
2009-11-20 21:50 - 2009-11-21 12:00 
2009-11-21 14:42 - 2009-11-23 00:50 
2009-12-29 02:21 - 2009-12-30 09:36 
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ity state recently analysed bv lBozzo et al.l ( l2010h . We re-analysed 
this same data-set, to properly compare it with our Suzaku obser- 
vation and check if the best fit found from Suzaku spectroscopy 
equally well describes the XMM — Newton spectrum, at a sim- 
ilar source luminosity state. IGR J08408~4503 was observed with 
XMM - Newton on 2007 May 29, with a live time of 35.7 ks 
(EPIC pn). 

Data were reprocessed using version 9.0 of the Science Analy- 
sis Software (SAS). Known hot, or flickering, pixels and electronic 
noise were rejected. Response and ancillary matrix files were gen- 
erated using the SAS tasks rmfgen and arfgen. We concentrate on 
the EPIC pn data: pn was operated in its Prime Full Window mode 
with all the CCDs in Imaging Mode. The EPIC pn observation used 
the thick filter. Source pn spectrum was extracted from a circular re- 
gion of 30" radius to avoid the edge of the CCD. We selected PAT- 
TERNS from to 4. Background counts were obtained from sim- 
ilar sized region offset from the source position. The background 
(selected with PATTERN=0 and above 10 keV) showed evidence 
of some flaring activity, which was filtered out, resulting in a net 
exposure time for the final source pn spectrum of 35.1 ks. The 
IGR J08408-4503 pn spectrum was not affected by pile-up, result- 
ing in an average net count rate of 0. 170±0.003 s~^. 

To ensure applicability of the x^ statistics, the net spectra 
were rebinned such that at least 30 counts per bin were present. 
All quoted uncertainties are given at 90% confidence level for one 
interesting parameter. 



3 ANALYSIS AND RESULTS 



Figure 2. Long-term Hght curve of the INTEGRAL Key Program data 
of the source region, with the 5cr upper-limits obtained with IBIS/ISGRI, 
overimposed on the times of the Suzaku observation. The IBIS/ISGRI up- 
per limits (18 — 60 keV) have been converted into flux in th e 1 — 100 keV 
range, adopting the average spectral parameters obtained bv lRomano et alj 
<2O09ll 



GRI (15 keV - 1 MeV; iLebrun et al.ll2003t) and PICsIT (175 keV- 
10 MeV: lLabanti et alj|2003h . An INTEGRAL pointing lasts about 
2ks. 

We analysed our IBIS/ISGRI private data between 55142.852 
MJD and 55195.400 MJD, where IGRJ08408-4503 was in the 
field of view (see T able [T|l, using the Of f-line Scientific Analysis 
package OS A 8.0 ( Goldwurm eTai]|2003h . 

IGR 108408^503 has never been detected, being always be- 
low the 5cr threshold of detection, both on timescales of a single ob- 
servation, each lasting 2 ks, and of a single satellite revolution (Rev. 
in Table 1). The source remains undetected even adding together 
the whole data-set of the seven revolutions analysed here (Table 1). 
In Figure |2] we show the 5a upper-limit lightcurve obtained with 
IBIS/ISGRI data, overimposed on the times of the Suzaku obser- 
vation, where we have converted the IBIS/ISGRI count rate, mea- 
sured in the 18 — 60 keV range, to the 1 — 100 ke V u nabsorbed flux, 
adopti ng the average spectral parameters obtained bv lRomano et al.l 
( l2009h . 

2.3 XMM - Newton 

We compare here our Suzaku results with a publicly available 
much shorter XAIM — Newton observation of a low luminos- 



3.1 Suzaku 

The background subtracted hght curve of IGR 108408^503 in two 
energy bands (softer and harder than 3 keV) together with their 
hardness ratio, are shown in Fig. [3] The hardness ratio does not 
show strong evidence of variability, thus we will study both the av- 
erage spectrum extracted from the whole exposure time, as well as 
time selected spectra from the initial part of the observation and 
from the peak of the two flares, separately. Fitting the whole hard- 
ness ratio as a function of time (with a bin time of 5 ks) to a constant 
model we obtain a value of 0.40 ± 0.03 (xl = 1.3 for 44 d.o.f.). 

A timing analysis was performed after having converted the 
event arrival times to the Solar System Barycentric frame. We 
searched for coherent periodicity, but found no evidence for pul- 
sations. We could place an upper limit on the pulsed fraction (de- 
fined as the semi-amplitude of the sinusoid al modulation devide d 
the mean count rate), computed according to lVaughan et al.l ( 119941) . 
of 40% at the 99% confidence level for periods in the range 16 s- 
500 s. 

We selected three spectra from the Suzaku data-set, as fol- 
lows: the first spectrum (which we call "persistent" spectrum in 
Table [2j was extracted accumulating all photons coming from the 
first 120 ks of the observation, with the lowest intensity level (see 
the light curve reported in Fig. |3); the second spectrum is from 
the peak of the first flare ("flare 1", between 145 ks and 160 ks in 
Fig.[3]l, while the third spectrum comes from the peak of the second 
flare ("flare 2", between 202 ks and 216 ks in Fig. [3). These time 
spans are subjective choices in order to extract the lowest intensity 
spectrum and the spectra at the peak of each flare. 

We first analysed the lowest intensity spectrum (the "persis- 
tent" pre-flare spectrum, with a net exposure time of 35 ks), always 
fitting together the spectra from the different XIS units (extracted 
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Figure 3. SuzakuDilSX net light curves of IGRJ084()8^503 in two energy ranges (the XISl count rate S is the background subtracted source rate in the 
range 0.4-3 keV, while H is the net count rate in the range 3-10 keV), together with the hardness ratio (H/S) versus time and versus the intensity (S+H). Bin 
time is 5 ks in the upper panels, while 10 ks in the lower panel. Start and stop times are in units of MJD - 40,000. 
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separately), with factors included in the spectral fitting to allow for 
normalization uncertainties between the instruments. A fit with an 
absorbed simple power law (PEGPWRLW in XSPEC) resulted in pos- 
itive residuals below 1 keV and in an unphysical very low energy 
absorption consistent with zero (resulting in a photon index of 1.5, 
with a reduced xi5/d.o.f.=l. 248/57). Thus in the following we will 
constrain the absorbing column density to be at least that predicted 
from optical m easurements towards the donor star (7Vh=3x10'^^ 
cm~'^: lLevder et al. 2007) . We then added a second component to 
the power law continuum (a thermal plasma model, MEKAL model 
in XSPEC), to account for the residuals evident at soft energies. 
We note that a similar low energy component was also found in 
the low intensity s tate observed in 2007 with XMM — Newton 
JBozzo et alj2010l) and explained by these authors as the likely X- 
ray emission from the blue supergiant companion. The double com- 
ponent model resulted in a much better fit (x^/d.o.f.=0.8 11/55), 
with the spectral parameters reported in Table |2] If we extrapo- 
late at high energies this best fit, we derive a flux of 4.35 x 10"^'' 
erg cm~^ s"'^ (15^0 keV), more than one order of magnitude 
lower than the upper limit we can place with Suzaku/HXD (see 
Sect. 12. It . If interpreted as X-ray emission from the companion 
star, this soft component should be obviously always present with 
similar spectral parameters, and with an X-ray luminosity con- 
sistent with that usually observed from OB supergiants (^10"^^- 
10^^ erg ~^: ICassinem et al.lll98ll) . Thus we included it in the fit- 
ting model also when analysing the other two spectra from the peak 
of the two flares. 

Fitting the first flare peak emission with this double compo- 
nent model (net exposure time of 7.8 ks; Table[2]l, we found a good 
fit (x^/d.o.f.= 1.041/71) but with a very low temperature fcTmckai 
<0.09 keV and a high normalization for this soft component, which 
translates into an X-ray luminosity (0.5-10 keV) of 4x 10"^^ erg ^^, 
too high to be compatible with the X-ray emission from the su- 
pergiant companion. Thus, we constrained the MEKAL tempera- 
ture and its normalization to be within the same variability ranges 
we obtained from the spectroscopy of the faint persistent spectrum 
(Table[2j, which we verified to be also consistent with the faintest 
emission observed in 2007 with XMM—Newton, and reanalysed 
here in Sect. 13.21 Constraining these two parameters in the double 
component continuum, we obtained a worse fit and wave-like resid- 
uals (x^/d.o.f.= 1.355/71) along the whole energy range. To solve 
this problem we added an additional absorption to the power law 
model, which can be explained with a different absorption towards 
the donor star and the compact object. This new double compo- 
nent model (PHABS X (MEKAL -l-PHABS X PEGPWRLW) adopting 
the XSPEC syntax) resulted in the best fit (xl/d.o.f.= 0.958/70) re- 
ported in Table |2] with an additional absorption, A^Hpow, towards 
the power law component of ~ 1.8x10^^ cm~^ and a power-law 
photon index, F, of ~2.1. The unabsorbed flux reported in Table |2] 
is the flux obtained correcting from the total absorption TVh, but 
not for A'^Hpow. Correcting also for A^'hpow, and measuring the un- 
absorbed flux coming only from the power law component, we ob- 
tained 6x10"^^ ergcm"^ s"^ (1-10 keV). 

Fitting the second flare peak emission with this same dou- 
ble component model (net exposure time of 2.7 ks; see Table |2] 
for the fit results), we obtained A'hpow^1.4x10^^ cm~^ for the 
power law component, consistent with that measured during the 
first flare. The unabsorbed flux measured from only the power law 
component is 7x10"'^'^ erg cm^^ s^^ (1-10 keV) which, extrap- 
olated to the range 15-40 keV, results into an unabsorbed flux of 
3xlO~^^ erg cm~'^ s~^, consistent with the upper limit placed 
from HXD observation. 




Energy (keV) 

Figure 4. Spectral fits to the Suzaku spectrum of IGRJ08408^503 ex- 
tracted from the whole observation. The upper panel shows the result adopt- 
ing an absorbed power law model (counts spectrum together with the wave- 
like residuals in units of standard deviation). The medium and lower panels 
display the best fit to the average IGR 108408^503 spectrum, obtained with 
a double component model composed by a MEKAL at low energies together 
with a differendy absorbed powerlaw model (PEGPWRLW in XSPEC). 



Finally, we performed a spectral fit to the average spectrum 
extracted from the whole Suzaku observation (net exposure time 
67.4 ks). The results obtained with the same double component 
model discussed before are reported in Table[2]and shown in Fig.|4] 
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Table 2. Results of the time resolved spectroscopy (persistent, flare 1 and flare 2) and of the average spectrum extracted from the whole SuzakulXlS 
observation. 



Spectrum 


A^H 


A^Hp 


ow 


r 


"'-'mekal 


normmekal 


Unabs Flux (1-10 keV) 


x2/d.o.f. 




(1022 cm-2) 


(1022 cm-2) 




(keV) 


(10-4) 


(10- 


-12 


erg cm 2 s i) 




Persistent 


o.mt'if 


_ 




1 t:7+0.18 
-■-^'-0.17 


22+°-°'' 


4 7+18.9 
^■'-2.7 






0.42 


0.811/55 




o-3o+°:r 


0.63 


+0.63 
-0.52 


r, |-,r, + 0.50 


23+''°'^ 
U-^-3_0.06 
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0.44 


0.752/54 
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- 
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<0.09 


100+iso 






9.2 


1.041/71 
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4.0 
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- 
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1.869/31 
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1.38 
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-0.75 
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0.29lHi 


23j:«, 
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1.492/30 


Average 


n fi9+0-20 


_ 




1 Qq+007 


0.18«:«2 


33111 






1.5 


1.646/265 




'^■^'-0.17 


1.28 


+0.24 
-0.23 


1-98-0.13 


0.23+ni 


11+- 






1.4 


1.211/264 



3.2 XMM - Newton 

The source low intensity state in the Suzaku observation is similar 
to that observed with XMM — Newton in 2007. Thus, we re- 
analysed here these same public EPIC pn data, in order to compare 
them with our Suzaku XIS results. In particular, we performed a 
different time selection spectroscopy with respect to that reported 
by Bozzo et al. (2010). We then reanalysed the spectra adopting 
the same best-fit model found from Suzaku analysis, to test if the 
same double component spectral model can explain both data-sets. 

The light curve observed in 2007 with XMM — Newton 
(EPIC pn data) is reported in Fig. |5] in two energy ranges (above 
and below 2 keV) together with their hardness ratio. In the same 
figure (lower panel) we also show an hardness-intensity plot, with 
a larger bin time, to clearly show that there is a clear correlation be- 
tween the hardness ratio and the source intensity, as already noted 
and discussed by Bozzo et al. (2010). 

Since hardness ratio and source intensity are correlated, an 
intensity-selected spectroscopy translates also into a hardness-ratio 
selected spectroscopy. In order to compare the Suzaku spectrum 
of the initial lowest intensity state ("persistent" emission in Table[2]l 
with the lowest intensity state observed with XMM — Newton, 
we extracted a spectrum where the EPIC pn rate (S+H in Fig|5] 
lower panel) was lower than 0. 14 s^^ (spectrum A in Table[3}. Dur- 
ing this low intensity, the contribution to the X-ray emission from 
the 0-type supergiant companion is maximum, thus allowing us to 
put better constraints to the companion soft energy component (e.g. 
ICassinelli et al.lll98lh . Three other intensity-selected spectra were 
extracted in the following EPIC pn total rate ranges (S+H in Fig|5] 
lower panel): 0.14 s-^-0.3 s"^ (spectrum B in Table[3](, 0.3 s^^- 
0.7 s"^ (spectrum C) and >0.7 s^^, to catch only the peaks of the 
flaring activity (spectrum D). 

We applied to the pn intensity-selected spectra the best fit 
model found from Suzaku data, a double component model with a 
continuum composed by a hot plasma model together with a power 
law model, differently absorbed (PHABS x (MEKAL +PHABS x 
PEGPWRLW) adopting the XSPEC syntax). We constrained the to- 
tal absorption (the first PHABS model, absorbing both continuum 
components) to be not less than the absorption towards the opti- 
cal companion (as during the Suzaku XIS analysis). We also con- 
strained the MEKAL parameters to be within the ranges found from 
the Suzaku XIS analysis. In Table [3] the spectral results of only 
the power law component with the total and the additional absorp- 
tion are reported. The fit results are good, demonstrating that the 
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nlime: 128.0 
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Bin time: 1024. 
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Figure 5. XMM - Newton EPIC pn light curves of IGR 108408^503 
in 2007, in two energy ranges (S is the total source rate in the range 0.2-2 
keV, while H is the total count rate in the range 2-12 keV), together with the 
hardness ratio (H/S) versus time and versus the intensity (S+H). Bin time 
is 128 s in the upper panels, while 1024 s in the lower panel. Start and stop 
times are in units of MJD - 40,000. 
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same model adopted to explain Suzaku spectra is a nice decon- 
volution of the XMM — Newton spectra too, although a simpler 
model (an absorbed power law) already results in a good descrip- 
tion of the spectra A and D. We note that the main responsible for 
the increasing of the hardness ratio during flares is very likely a 
harder power law high energy component rather than an increasing 
of the absorbing column density (which remains constant, within 
the uncertainties), as already reported by Bozzo et al. (2010), but 
adopting a different deconvolution of the spectrum. 



In the framew ork of the clumpy wind model developed by 
iDucci et alj ( 120091) , we can derive the distance of the neutron star 
from the supergiant companion, from the durations and luminosi- 
ties of the flares observed with Suzaku, a nd adopting the ex pan- 
sion law for the clump s ize (see Eg. [9] in Ducci et alj|2009l ; see 
also iRomano et al.ll2009l and ISidoli et alj|2009l for previous appli- 
cations of this method to IGRJ08408-4503). Assuming a wind 
velocity of v^ — 1800 km s~^ and a neutron star with mass 
Mx — 1.4 Mq, the accretion radius is: 



Ra. 



2GM, 







1.2 X lO"^ cm 



(1) 



4 DISCUSSION AND CONCLUSIONS 

We report here on the longest observation of a low luminosity state 
of a SFXT, spanning about three consecutive days (with the ob- 
vious gaps due to the satellite orbits), resulting in a net exposure 
time of 67.4 ks. The Suzaku IGR J08408^503 observation un- 
veils a source luminosity level at about 4x10^^ erg s^^ for the 
first ~120 ks, followed by two long flares (with a duration of about 
45 ks each) with a count rate a factor of ^--^3 higher than the initial 
"persistent" pre-flare emission. The source has not been convinc- 
ingly detected with HXD and we could place the best upper limit 
to date to the hard X-ray emission (15-40 keV) of a SFXT in the 
low intensity state, at a level of 6x 10^^ erg s~^. 

The time selected spectroscopy of the XIS observation re- 
sulted in a best fit model composed by a double component con- 
tinuum (a soft hot plasma model together with a power law domi- 
nating at high energies), differently absorbed, with the presence of 
an extra absorption of the high energy power law component. This 
additional absorbing matter, of the order of ~l-2x lO'^^ cm~^, is 
very likely local to the compact object, since it absorbs only the 
power law component. It could be due to dense wind clumps pass- 
ing in front of the compact object and close to it, or more likely to 
the accreting matter itself, since this additional absorption is more 
clearly detected during the two long duration flares. The spectral 
parameters of the double component model suggest that the low 
energy soft component is produced by the supergiant companion, 
while the high energy component is too hard for a stellar origin (we 
will not repeat here the convincing arguments already discussed by 
Bozzo et al. 2010). Both the temporal and spectral properties indi- 
cate that in IGRJ08408-4503 the compact object is still accreting 
matter even at the low intensity state we observed with Suzaku: in- 
deed, also in the lowest luminosity state of 4x 10"'^ erg s~^, during 
the first part of the observation (long term pre-flare emission), the 
Suzaku spectrum is hard and power-law-like at energies higher 
than 1 keV, with a photon index, F, of 1-2. Similar photon in- 
dexes have been measured in the 1-10 keV spectra of four other 
SFXT s monitored with S wift during their long-term intermediate 
states JSidoli et al]|2008l) , although at higher luminosities (10^^- 
10^* erg s~^). The presence of flares, together with the hard X- 
ray emission, is strongly indicative of the fact that accretion on the 
compact object is still ongoing in IGR 108408-4503, even at very 
low luminosities. 

SFXTs behaviour is still not explained (both the outbursts 
and the high dynamic range), although several possible mecha- 
nisms have bee n proposed since th e times of their discovery with 
INTEGRAL ('Sguera et al.'2005'). Among the different proposed 
explanations (see Sidoli 2009 for a review) which involve the wind 
properties and/or the characteristics of the compact object (neutron 
star magnetic field and spin period) there is a general consensus on 
the fact th at the supergian t wind is very likely not homogeneous 
("clumpy" : lin't Zanj|2005h . 



From the time duration of the flares, ta ~ 4.5 x 10 s, we obtain 
the radius of the clumps accreted by the neutron star: 



Rci = 



'. 4 x 10^^ cm 



(2) 



Extrapolating flux of the two flares to a broader energy range, 
Fx ~ 1 X 10^^^ erg cm^^ s^^ (1 — 40 keV), and assuming a 
distance of d = 3 kpc, we obtain a flare X— ray luminosity L^ ~ 
10'''* erg s~*. Hence, from the accretion X-ray luminosity and the 
continuity equation A/accr ~ pv^^nRii, we obtain a clump density 
of Pel ~ 6.6 X 10"*'' g cm~^, and a mass of Md ^ 10^^ g. 
For the supergiant companion we assumed the following physical 
properties: a radius Rob = 23.8 Rq, a mass AfoB = 30 M0, a 
luminosit y Ior L/Lp, = 5.8 47 and an effective temperature T^s ~ 
34000 K jVaccaetal.ll996l) . We thus obtain, from the clump mass- 
radius equation of Ducci et al. 2009 (Eq. [20]), the starting radius of 
the clump: Rci,i — 1.7 x 10** cm. Assuming that the sonic radius 
Rs equals the supergiant radius, and a starting wind velocity at Rs 
of Uw(-Rs) ~ 3 x 10*^ cm s~* (Bouret et al. 2005), we obtain, from 
the expansion law of the clump, a distance of the neutron star of 
r ~ 1 0*^ cm, consistent with our previous findings JRomano et al.l 
I2OO9I) . 

The orbital separation derived here implies an orbital period 
of about 35 days, which is somehow intermediate between the nar- 
row orbits of o ther SFXTs (e.g., 3.3 days in IGR 116479-4514, 
Ijain et alj2009n a nd the lon gest orbit al period measeured to date in 
a SFXT (165 davs. lSidoU eT al. 2006, Romano et al.l2009l) . whereas 
it is very si milar to tha t measured in SFXTs like, e.g, SAX J18 18.6- 
1703 (Zurita Heras & Chatv 2009, Bird et al. 2009) . Even the long 
duration low luminosity emission at a few lO^'^ erg s~* (c lose to the 
quiescent emission observed, e.g., in IGR J17544-2619. lin't Zand 
2005) is consistent with a wide orbit, in order to allow accretion 
at a low level. Also the long term light curve reported here from 
INTEGRAL observations seems to point to a long duration of 
the low intensity emission in IGR 108408^503. 

The presence of short flaring activity as observed with 
XMM - Newton in IGR 108408^503 during a low intensity 
state has already been caught in oth er SFXTs (e.g., in IGR 117544- 
2619. iGonzalez-Riestra et alj|2004l . their Fig. 3) with flares dura- 
tions of a few hundred seconds. On the other hand, the Suzaku 
pointing spanning three days allowed us for the first time to probe 
the source light curve on much longer timescales, catching longer 
flares (~45 ks) than usual, with peak rates which are only ~3 times 
higher than the average pre-flare intensity. The observed proper- 
ties of these flares, their long duration together with the low dy- 
namic range with respect to the pre-flare emission, suggested us 
to search for other possible explanations, although in the frame- 
work of the structure of the supergiant wind. A viable alterna- 
tive to the accretion of single clumps could be that the flares are 
produced from the accretion of large scale wind structures, like 
gas streams, or the so-called corotating interaction regions (CIRs, 
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Table 3. XMM — Newton results of the source intensity (and hardness) selected spectroscopy (EPIC pn) adopting the model PHABS X (MEKAL +PHABS 
X PEGPWRLW). See Sect.|3]2]for details. 



Spectrum A^h A^Hpow 

(10^2 cm-2) (io22 cm-2) 



Observed Flux (1-10 keV) X^/d.o.f. 
(10~^3 ergcm~2 g^i) 



A 


0.67; 


B 


0.61; 


C 


0.49; 


D 


0.50 



+0.06 
-0.37 
+0.13 
-0.15 
+0.33 
-0.18 
+0.10 
-0.10 



<2.1 

,+0.38 
-0.32 

"si+^-so 



o.sr 



2.94 
2.42 
1.95 
1.46 



+0.92 
0.69 

+0.21 
0.19 

+0.25 
0.23 

+0.25 
0.27 



1.8 


1.209/20 


6.4 


1.090/64 


23 


0.993/19 


58 


0.713/11 



lMullanlll984n . CIRs form in a stellar wind when the rotating star 
emits wind in a non-spherically symmetric manne r. These wind 
structures were originally studied in the solar corona ( IMullanl 19840 
and form as interaction regions between slow and fast stream- 
lines. They were later suggested to be present also in rotating O- 
star winds iCranmer & Owockilll99 6). with a mild density con- 
trast of a factor of 2-3 with respect to the undisturbed wind. These 
structures hav e a typical thickness of ~0. Ir at radial distance r 
( lMullaall984h . Thus, assuming a radial distance r ~ 10^^ cm 
and a stream thickness, AS~0. Ir, we can derive the neutron star 
velocity, w„s from AS and the flare duration At/iarc (AS=v„s x 
At/jare)- Assumiug Atfiare='^5 ks and AS~0.1r-=10^^ cm, we 
find v„s=200 km s~^, which is consistent with the neutron star or- 
bital velocity in a 35 days orbit around a 30 Mq companion star. 



ACKNOWLEDGMENTS 

This work was supported in Italy by contract ASI/INAF 1/088/06/0. 
P.E. acknowledges financial support from the Autonomous Re- 
gion of Sardinia through a research grant under the program PO 
Sardegna FSE 2007-2013, L.R. 7/2007 "Promoting scientific re- 
search and innovation technology in Sardinia". This work is based 
on data from observations with Suzaku, INTEGRAL, and 
XMM — Newton. Suzaku is a Japan's mission developed at the 
Institute of Space and Astronautical Science of Japan Aerospace 
Exploration Agency in collaboration with U.S. (NASA/GSFC, 
MIT) and Japanese institutions. XMM — Newton is an ESA sci- 
ence mission with instruments and contributions directly funded 
by ESA Member States and the USA (NASA). INTEGRAL is 
an ESA project with instruments and science data centre funded 
by ESA member states (especially the PI countries: Denmark, 
France, Germany, Italy, Switzerland, Spain), Czech Republic and 
Poland, and with the participation of Russia and the USA. This re- 
search has made use of HEASARC online services, supported by 
NASA/GSFC. 



REFERENCES 

Barba R., Gamen R., Morrell N., 2006, Astron. Tel., 819, 1 
Bird A. J., Bazzano A., Hill A. B., McBride V. A., Sguera V., 

Shaw S. E., Watkins H. J., 2009, MNRAS, 393, LI 1 
Bouret J.-C, Lanz T, HilherD. J., 2005, A&A, 438, 301 
Bozzo E., Stella L., Ferrigno C, Giunta A., Falanga M., Campana 

S., Israel G., Leyder J. C, 2010, ArXiv 1004.2059 
Cassinelli J. P, Waldron W. L., Sanders W. T, Harnden Jr. F. R., 

Rosner R., Vaiana G. S., 1981, ApJ, 250, 677 
Cranmer S. R., Owocki S. P, 1996, ApJ, 462, 469 



Ducci L., Sidoh L., Mereghetti S., Paizis A., Romano P., 2009, 

MNRAS, 398, 2152 
Fukazawa Y., Mizuno T, Watanabe S., Kokubun M., Takahashi 

H., Kawano N., Nishino S., Sasada M., et al. 2009, PASJ, 61, 17 
Goldwurm A., David P., Foschini L., Gros A., Laurent P., 

Sauvageon A., Bird A. J., Lerusse L., et al. 2003, A&A, 411, 

L223 
Gonzalez-Riestra R., Oosterbroek T, Kuulkers E., Orr A., Parmar 

A. N., 2004, A&A, 420, 589 
Gotz D., Falanga M., Senziani F, De Luca A., Schanne S., von 

Kienlin A., 2007, ApJL, 655, LlOl 
Gotz D., Schanne S., Rodriguez J., Leyder J.-C, von Kienlin A., 

Mowlavi N., Mereghetti S., 2006, Astron. Tel., 813, 1 
Gruber D. E., Matteson J. L., Peterson L. E., Jung G. V., 1999, 

ApJ, 520, 124 
in't Zand J. J. M., 2005, A&A, 441, LI 
Jain C, Paul B., Dutta A., 2009, MNRAS, 397, LI I 
Kennea J. A., Campana S., 2006, Astron. Tel., 818, 1 
Koyama K., Tsunemi H., Dotani T, Bautz M. W., Hayashida K., 

Tsuru T G., Matsumoto H., Ogawara Y, et al. 2007, PASJ, 59, 

23 
Labanti C, Di Cocco G., Ferro G., Gianotti F, Mauri A., Rossi 

E., Stephen J. B., Traci A., Trifoglio M., 2003, A&A, 411, L149 
Lebrun F, Leray J. P., Lavocat P., Cretolle J., Arques M., Blondel 

C, Bonnin C, Bouere A., et al. 2003, A&A, 411, L141 
Leyder J.-C, Walter R., Lazos M., Masetti N., Produit N., 2007, 

A&A, 465, L35 
Lu F J., Aschenbach B., 2000, A&A, 362, 1083 
Lund N., Budtz-j0rgensen C, Westergaard N. J., Brandt S., Ras- 

mussen I. L., Hornstrup A., et all 2003, A&A, 411, L231 
Masetti N., Bassani L., Bazzano A., Dean A. J., Stephen J. B., 

Walter R., 2006, Astron. Tel., 815, 1 
Mereghetti S., Sidoh L., Paizis A., Gotz D., 2006, Astron. Tel., 

814, 1 
Mitsuda K., Bautz M., Inoue H., Kelley R. L., Koyama K., Ku- 

nieda H., Makishima K., Ogawara Y, et al. 2007, PASJ, 59, 1 
Moretti A., Pagani C, Cusumano G., Campana S., Perri M., 

Abbey A., Ajello M., Beardmore A. P, et al. 2009, A&A, 493, 

501 
Mullan D. J., 1984, ApJ, 283, 303 
Romano P., Sidoli L., Cusumano G., Evans P. A., Ducci L., Krimm 

H. A., Vercellone S., Page K. L., et al. 2009, MNRAS, 392, 45 
Romano P., Sidoli L., Cusumano G., La Parola V., Vercellone S., 

Pagani C, Ducci L., Mangano V., Cummings J., Krimm H. A., 

Guidorzi C, Kennea J. A., Hoversten E. A., Burrows D. N., 

Gehrels N., 2009, MNRAS, 399, 2021 
Romano P., Sidoli L., Cusumano G., Vercellone S., Mangano V, 

Krimm H. A., 2009, ApJ, 696, 2068 



Suzaku observes IGR J08408-4503 during a long low luminosity state 9 



Romano P., Sidoli L., Mangano V., Mereghetti S., Cusumano G., 

2007, A&A, 469, L5 
Romano R, Sidoli L., Mangano V., Vercellone S., Kennea J. A., 

Cusumano G., Krimm H. A., Burrows D. N., Gehrels N., 2008, 

ApJL, 680, LI 37 
Sguera V., Barlow E. J., Bird A. J., Clark D. J., Dean A. J., Hill 

A. B., Moran L., Shaw S. E., Willis D. R., Bazzano A., Ubertini 

R, Malizia A., 2005, A&A, 444, 221 
Sguera V., Bazzano A., Bird A. J., Dean A. J., Ubertini R, Barlow 

E. J., Bassani L., Clark D. J., Hill A. B., Malizia A., Molina M., 

Stephen J. B., 2006, ApJ, 646, 452 
Sidoli L., 2009, Advances in Space Research, 43, 1464 
Sidoli L., Paizis A., Mereghetti S., 2006, A&A, 450, L9 
Sidoli L., Romano P., Ducci L., Paizis A., Cusumano G., Mangano 

v., Krimm H. A., Vercellone S., Burrows D. N., Kennea J. A., 

Gehrels N., 2009, MNRAS, 397, 1528 
Sidoli L., Romano P., Mangano V., Cusumano G., Vercellone S., 

Kennea J. A., Paizis A., Krimm H. A., Burrows D. N., Gehrels 

N., 2009, ApJ, 690, 120 
Sidoli L., Romano P., Mangano V, Pellizzoni A., Kennea J. A., 

Cusumano G., Vercellone S., Paizis A., Burrows D. N., Gehrels 

N., 2008, ApJ, 687, 1230 
Smith D. M., Heindl W. A., Markwardt C. B., Swank J. H., 

Negueruela I., Harrison T. E., Huss L., 2006, ApJ, 638, 974 
Takahashi T, Abe K., Endo M., Endo Y., Ezoe Y., Fukazawa Y, 

Hamaya M., Hirakuri S., et al. 2007, PASJ, 59, 35 
Ubertini P., Lebrun R, Di Cocco G., Bazzano A., Bird A. J., 

Broenstad K., Goldwurm A., La Rosa G., et al. 2003, A&A, 411, 

L131 
Vacca W. D., Garmany C. D., Shull J. M., 1996, ApJ, 460, 914 
Vaughan B. A., van der Klis M., Wood K. S., Norris J. P., Hertz 

P., Michelson P. R, van Paradijs J., Lewin W. H. G., Mitsuda K., 

Penninx W., 1994, ApJ, 435, 362 
Vedrenne G., Roques J.-R, Schonfelder V, Mandrou P., Lichti 

G. G., von Kienlin A., Cordier B., Schanne S., et al. 2003, A&A, 

411, L63 
Winkler C, Courvoisier T., Di Cocco G., Gehrels N., Gimenez A., 

Grebenev S., Hermsen W., Mas-Hesse J. M., et al. 2003, A&A, 

411, LI 
Zurita Heras J. A., Chaty S., 2009, A&A, 493, LI 

This paper has been typeset from a Tj5£/ WT^Si file prepared by the 
author. 



